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Abstract 

State tomography on qubit pairs is routinely carried out by measuring the two qubits separately, while one 
expects a higher efficiency from tomography with highly symmetric joint measurements of both qubits. Our 
numerical study of simulated experiments does not support such expectations. 
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1. Introduction 

Quantum state tomography makes efficient use of the resources if it employs a probability operator 
measurement (POM) with the minimal number of outcomes — 4 = 2 2 for single-qubit tomography, 16 = 4 2 
for qubit-pair tomography, d 2 for tomography in a d-dimensional Hilbert space. Ideally, one would like to 
acquire tomographic data with the aid of a symmetric informationally complete (SIC) measurement P, 0] 
because it probes the state space in an unbiased way and promises minimal statistical errors on average 
when many, uniformly distributed, quantum states are measured [H, HI- 

The SIC tetrahedron measurement 0, [B[ has been implemented for polarization qubits of single photons 
and the product POM that is composed of the joint measurements of two single-qubit SIC POMs has 
been used successfully for two-qubit state tomography This product POM is, however, not a SIC POM 
for the qubit pair, which invites the question whether the realization of the two-qubit SIC POM, with its 
challenge of implementing projections to entangled states, is worth the effort. Our answer is: Hardly. 

The evidence that justifies this answer is numerical. We generate simulated tomographic data, both for 
the product POM and the SIC POM, for two-qubit states that are chosen either specifically or randomly. 
Then we estimate the state from the simulated data and compare the reconstructed state with the pre-chosen 
true state. The figure of merit is the trace-class distance between the two states. 

The paper is organized as follows. We introduce notation and terminology in Sec. [2] where we briefly 
review schemes for state tomography. Then, in Sec. El we remark on group covariant SIC POMs and specify 
the product POM and SIC POM that we use for the numerical study. Its ingredients — state estimation 
and the generation of random samples of states — are described in Sees. [Hand respectively. The results 
of the simulated experiments are presented in Sec. E]for tomography of maximally entangled states, and in 
Sec. [7] for averages over random samples of states. We close with offering our conclusion in Sec. [H 
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2. Measurement schemes for state tomography 

is a decomposition of the identity into J probability operators IX,- , which we call outcomes, 

^2 n j ■ = 1 with n j ■ > ■ (1) 

3=0 

The probabilities 

.7-1 

Pj = tr {pllj} with /Jft = 1 (2) 

3=0 

permit a reconstruction of the statistical operator p, the state of the physical system, if the POM is informa- 
tionally complete (IC). In a d-dimensional Hilbert space, an IC POM must have at least d 2 outcomes, and 
a minimal IC POM has exactly J = d 2 outcomes and no redundant ones, and then the state p is uniquely 
expressed in terms of the probabilities pj and the reconstruction operators Pj that are dual to the outcomes 
II.-- 

d 2 -i 

P=Y1 p i P 3 ' where tr{P.,n fc } = 6 jk ■ (3) 
j=o 

A SIC POM is a minimal IC POM with the symmetry property that all outcomes are unitarily equivalent 
with "equal angles" between each pair of outcomes, 

tr{n.,n fc } = aSjk + b with tr{ITj} = a + bd 2 = ^ 

and l<tr{n j 2 }= a + 6<l. (4) 
The reconstruction operators for such a SIC POM are 

Pi = i(lL,--M), (5) 

so that 

d 2 -l 

tr{P J } = l, tr{Pj-Pfc} = -((Sj-fc - bd) , Y, P i= d - ( 6 ) 

3=0 

These translate d _1 < tr{p 2 } < 1 into 

d 2 -l 
3=0 

a second constraint obeyed by the probabilities pj in addition to the unit sum in {5]) . The equal sign in Q 
applies when p is a pure state, p 2 = p, and only then. 

When the d 2 outcomes 11^ of a SIC POM are rank-1 operators, that is: when they are subnormalized 
projectors to pure states, then a + b assumes the maximal value permitted by (J3J), 

rank-1 SIC POM: a + 6= -L, a = bd = ' 



d 2 ' d + d 2 ' 

n 2 = i nj , P j = (d + d 2 )U j -l. (8) 

Almost all SIC POMs studied in the published literature are of this rank-1 kind. One exception is the 
two-qubit SIC POM of Ref. whose rank-3 outcomes are optimal entanglement witnesses. In the present 
context, we will only deal with standard rank-1 SIC POMs. 



1 POM, with its emphasis on "probability" and "measurement" is quantum physics jargon. The corresponding popular 
mathematical term POVM (positive-operator valued measure) makes reference to measure theory. 
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3. Group covariant SIC POMs and their fiducial states 

As conjectured in [H, Q, for every Hilbert space dimension d, there exist rank-1 SIC POMs which are 
covariant with respect to the Heisenberg-Weyl group, also called the generalized Pauli group because it is 
the G?-dimensional analog of the d = 2 group of unitary transformations effected by the Pauli operators a x , 
o z and their products, including a x a z — —i<J y as well as the identity a 2 = 1. Indeed, SIC POMs of this 
group covariant kind are typical. 

The Heisenberg-Weyl group refers to a pre-chosen basis of kets |0), |1), . . . , \d— 1). They are eigenkets 
of the d-periodic unitary operator Z, the analog of o z , 

Z\k) = |fc)e i27rfe / d , Z d = l, (9) 

and are cyclically permuted by the rf-periodic unitary operator X, the analog of a x , 

X\k) = \k + 1) forfc = 0,l,...,d-2, X|d— 1) = |0) , X d = 1 . (10) 

The (m, n)th element in the abelian Heisenberg-Weyl group of unitary transformations is the mapping 

^ X m z n F Z~ n X~ rn Z n X m F X~ m Z~ n 

for m,n = 0, 1, . . . , d — 1 , (11) 

where F is any linear operator on the d-dimensional Hilbert space of kets. Although X and Z do not 
commute, ZX = the order of factors does not matter in the transformation (fTT|) . and so we have 

d 2 group elements. 

A group covariant SIC POM is composed of the outcomes 

n m „ = x m z n u 0Q z- n x- m , (12) 

where the d 2 pairs m, n take over the role of label j in Q and Hoo is the fiducial outcome. The ergodic 
property of the Heisenberg-Weyl group, 

d-l 

X m Z n FZ- n X- m =dtr{F}, (13) 

m,n— 

ensures that the n mn s are proper outcomes of a POM if Hoo is a probability operator with tr-{Tloo} = 1/d. 
But the rank-1 SIC POM condition of (gj with a and b from JHJ) is obeyed only if n 00 = |f)d _1 (f| is such 
that 

|(f|A m Z"|f)| 2 = 1+d ^Jn forTO;71 = 0;lj ... 5d _ lj (14) 

where the normalized |f) is the fiducial ket or seed The set of equations (fl4|) can serve as the basis for 
a numerical search, successfully completed by Renes et al. up to d = 45 [2| , but there are also analytically 
known fiducial kets and their rank-1 POMs. 

For the qubit case d = 2, we have the familiar tetrahedron measurement with the outcomes 



1 



To = 


4 v / 3 




1 


Ti = 


4 v / 3 




1 


T 2 = 


4v/3 




1 


T 3 = 


4\/3 



3 + a x + o-y + a z ) , 

3 + o x - o- y - er 2 ) = o- x T a a x , 

3 - o x - o-y + cr 2 ) = a z T Q a z , 

3 - o x + o-y - er 2 ) = o- x cr z T a a z a x , (15) 
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for which 



E T i = 1 ' tr {^} = ^ tr{T i ^} = i ( 5 ifc + i. 



Its name derives from the geometry of the four Pauli vectors tj , 



T 4 = -(1 + U-a) 



-_4 1 
' ~~ 3 _ 3 ' 



(16) 



(17) 



whose tips are the corners of a tetrahedron inside the Bloch sphere. See Ref. p>] for the properties of 
the tetrahedron measurement, and Ref. Q for its implementation for experimental tomography of photon 
polarization qubits. 

The tensor product of two tetrahedron SIC POMs, 



T„ 



(18) 



m,n— 



is a perfectly suitable POM for two-qubit tomography. It is central to the Singapore protocol for quantum 
key distribution [ltj| . and has been used in the experiments by Ling et al. Q. This product POM is, however, 
not a SIC POM for the qubit pair. 

For two-qubit SIC tomography, we can use any one of the known SIC POMs for d = 4, but they are not 
equivalent in their two-qubit properties, such as the concurrence of the outcomes llj . We prefer a particular 
SIC POM, generated from Appelby's fiducial state 



|f A ) = (|00),|01),|10),|ll)) 



2V3TG 



( l + e~ i7r /4 \ 

e i7r/4 + iG _3/2 

1 - e~ i7r / 4 

y e i7r/4_ iG -3/2 J 



(19) 



where G — (V5 — l)/2 is the golden ratio. The transformations of the Heisenberg-Weyl group refer to the 
two-qubit kets \j1j2) = \ji) <8> lj'2) in the order stated, so that we have 



Z = — ^ a z (g) (1 - icr z ) , X = -(1 + a x ) <S> a x - ^(1 - a x 

for the two-qubit versions of the basic unitary operators in ^ and (fTTTj) . 

The outcomes II m „ of the particular SIC POM that we get from the seed {T 



(20) 



n, 



x m z n u A )-tf A \z- n x- 

d 



for m, n = 0, 1, 2, 3 , 



(21) 



possess an additional symmetry: all the 16 two-qubit states dll mn have the same concurrence of -y/2/5. 
These two-qubit states are typical in the sense that their squared concurrence equals the average squared 
concurrence of all pure two-qubit states [l3T | . As a consequence of the common concurrence, one can turn 
the II mn s into each other with the aid of local unitary transformations that act on one of the qubits only. 
This is a rather peculiar feature, not shared by other SIC POMs for qubit pairs 11 1. 



4. Techniques for state estimation 

The data acquired in the simulated experiments are simply the counts of clicks of each of the 16 detectors, 
one for each outcome of the two-qubit POMs — the product POM of (JTBJ) and the SIC POM of |TS| and (f2"T|) . 
Owing to the statistical fluctuations that originate in the quantum indeterminism, the relative frequencies 
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of detector counts are not equal to the probabilities Pm'n^ = tr{pT mn } and pm^ — tr{pll m „}, at best the 
respective relative frequencies pm™ d ' and pmn approximate the intrinsic probabilities. 
Therefore, there is danger in relying on (|3|). The reconstructed states 



P c r d) = E P { £n d) (ST m - 1) <g> (6T n - 1) , 

m,n— 

Pcst IC) = E pL S n C) (20n„ m -l) (22) 

m,n— 

are not guaranteed to be positive operators. It is possible, and indeed happens regularly when the true 
state has high purity, that these raw data (RD) estimates have negative eigenvalues because the relative 
frequencies do not obey restrictions such as the appropriate version of ([7]). In the context of singlc-qubit 
tomography with the tetrahedron measurement, this matter is discussed in section IV C of Ref. A 
thorough cov erag e of many aspects of quantum state estimation is given in the 2004 monograph of Paris 
and Rehacek [1J|. 

The probability that a total of N detector clicks is recorded with the actual relative frequencies p mn is 
the likelihood C(p). Its logarithm is 

3 

log(£(p)) =N E 

p mn logp mn , (23) 

m,n— 

where the unknown intrinsic probabilities p m n derive from the true statistical operator p that we wish to 
estimate. In the maximum likelihood (ML) strategy we take the p for which C(p) is largest as our best guess 
for the true p, 

C(p C st) = max£(p) . (24) 
p 

The ML estimate is identical with the respective RD estimate of (|22p if the RD estimate is permissible, that 
is: if there is a physical p whose probabilities are equal to the recorded relative frequencies. Accordingly, 
we need to calculate the ML estimate only if the RD reconstruction (|2"2"|) fails. 

The ML estimate is found by a numerical iteration that begins with an arbitrary guess p'fsl, for which 
the completely mixed state p^ s l — 1/4 is an option, and then determines improved estimates successively in 
accordance with the update rule 

(fe+i) _ [l + e k Rk]pjsl [l + ffc-Rfc] , , 

tr[[l + e k R k ]pi:l[l + e k R k ]\ 

where (replace IT m „ by T mn for the product POM) 

3 - „ 

R *= E ; r (fc)n x - 1 ( 2fi ) 

m,n— ^r< p cs t limn f 

and efc > is a suitably chosen increment (see below). The iteration stops when the trace class distance 
trj \R k Pesl 1 1 is below the pre-set accuracy threshold. The rule ([23]) . which is a variant of the recipe given 

in chapter 3 of [3], is such that all Pest s are assuredly positive and normalized to unit trace. 
In a bit of detail, the important steps of the iteration are as follows. 

Given p^ st ; , 

i. Compute R k ; terminate the iteration if tr| \R k p^ t 1 1 is below the accuracy threshold, oth- 
erwise proceed. 
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ii. Use two trial values for e^- to compute two p[st s an d determine the value of the likelihood 
£(^ t +1) ) for both. 

iii. Combine these two with the e& = value £(p^ t ') and compute a quadratic function of e& 
that interpolates between the three support values. 

iv. Find the value for which the quadratic function assumes its maximum. 

v. Use this maximizing ej, in (|25p for the update p^l — > p^ 1 ^. 
vi. Repeat. 

With the optimization of Ck in steps ii-iv, convergence is quite fast in practice. 

5. Choosing states at random 

For the simulated state tomography, we choose the true states either specifically or at random. We 



generate the random states by a procedure described in Refs. jl3l. 1 151]. A random rank-r state p is selected 
by first populating the entries of a 4 x r matrix Y such that the real and imaginary parts follow noncentral 
Gaussian distributions offset by mean values given by the matrix M = E(Y). Then 



p = (|00>, |01>, |10>, lll))-^^- 



/ <oo| \ 

(01| 
(10| 

1(111/ 



(27) 



is a rank-r state in the random sample that is characterized by matrix M. The choice M — for r = 4 
gives a uniform sample with respect to the Hilbert-Schmidt distance, by other choices we can regulate the 
average purity of the states in the random sample. 

Random samples of pure states (r — 1) are unitarily invariant and thus uniform. So are random samples 
of unbiased mixed states (M — 0), but if we choose M ^ for a mixed states sample, there will be a bias 
because we lose the unitary equivalence. The average performance of the SIC POM will depend on the 
sample or, equivalently, which one of the unitarily equivalent SIC POMs we use for the same biased sample. 

Once the true state p is chosen for the simulation, we use the resulting quantum probabilities in con- 
junction with a standard (pseudo-)random number generator to get a total of N simulated detector clicks. 
The relative frequencies p mn of detector clicks then determine the RD estimate of p and, if need be, also 
the ML estimate. The quality of the estimate is judged by the trace-class distance, 

D = hv{\p ~ Pcst \} , (28) 

with < D < 1. We have D = only if the two states are equal, p cs t — p, and the maximal distance D = 1 
when they are orthogonal p es tP = 0. Typically, the distance decreases D cx N^ 1 / 2 as the sample size N 
increases. 



6. Tomography of maximally entangled pure states 

We now compare the tomographic performances of RD and ML reconstruction for pure states for which 
the failure of RD estimation is typical rather than exceptional. In particular, we choose as the true states 
the four standard Bell states that are given by 
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Figure 1: Simulated tomography of two-qubit Bell states with the product POM (crosses + and solid line) or the SIC POM 
(crosses X and dotted line) with raw-data state estimation (left-side plots) and maximum-likelihood state estimation (right-side 
plots). The trace-class distance D aV g is averaged over a hundred simulation runs with up to 6000 measured qubit pairs, all 
prepared in the respective Bell state. The solid and dotted lines are least-square fits to a power law D a vg = a/N c . 
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I^H-LdOl^llO)) and |$ ± ) = i=(|00)±|ll)). (30) 

The plots in the left column of Fig. [1] show the N dependence of the distance between the RD estimate 
and the respective Bell state. Not unexpectedly, the SIC POM outperforms the product POM. But this 
observation is, in fact, misleading as almost all RD estimates are unphysical. The comparison with the ML 
estimates in the right column is more to the point, and there the picture is reversed: rather unexpectedly, 
the product POM gives better estimates than the SIC POM. 

The simulated product-POM data compare well with the experimental data of Ref. 0], where a RD 
estimate was used. The experimental values for the average trace-class distance D a vg showed a dependence 
cx TV -1 / 2 on the total number of detector clicks, and our simulated data follow the same reciprocal square- 
root law. 

The striking similarity among the four plots of each column in Fig. [T] is not accidental but could in fact 
have been anticipated. The transformations effected by the squares of the Heisenberg-Weyl operators X 
and Z of (f20f interchange the Bell states, 

X 2 p Q X 2 =p x , Z 2 p x Z 2 = p y , X 2 Py X 2 = p z , Z 2 p z Z 2 = p , (31) 

where X 2 = a x ® 1 and Z 2 = 1 ® a z act on one of the qubits only. It follows that the respective sets of 
probabilities p mn of the Bell states are permutations of each other. 

For a more systematic quantitative comparison between the product POM and the SIC POM, we intro- 
duce a performance factor rj, guided by the usual understanding of optimal tomography with SIC POMs. 
Suppose the power laws have the form D = a/N c , and jV( prod ) and iV*( SIC ) are the respective number of 
qubit pairs that need to be measured to reach the desired threshold value Z?thr- Then 

_ ^od) _ /^od) X l/e<- d > f ^ ^ 1/<P*> 

77 ~ iV(sic) - ^ D thr ) U (SIC V 

defines r\. For r\ > 1, the larger its value, the faster the convergence rate of the SIC POM power law 
compared to that of the product POM for a fixed benchmark distance D t hr- Similarly, ij < 1 indicates that 
the product POM is more efficient. 

From top to bottom, the r\ values for the data in Fig.[T]are 1.32, 1.31, 1.37, and 1.29 for the RD estimates 
of the left column, and 0.42, 0.43, 0.41, and 0.42 for the ML estimates of the right column, all referring to the 
benchmark value -Dthr = 0.1. It is staggering how much better the product POM performs as soon as proper 
ML estimation of the state is incorporated. This observation is not limited to the maximally entangled Bell 
states. We have seen it also for quite a few randomly chosen pure two-qubit states and for the three-qubit 
W state |W) = ( 1 001) + |010) + 1 100)) /-s/3- A more extensive study of three-quit tomography appears to be 
worthwhile. 

That said, we must also note that the four Bell states of (f2"9"|) are not typical among the maximally en- 
tangled states, owing to their particular alignment with the outcomes of the two POMs under consideration. 
When one averages over a large sample of maximally entangled states, see the bottom right subtable in 
Table [TJ the picture is different: the average r\ value is a bit larger than 1. 



7. Averaged performances 

The data for average performances of the product POM and the SIC POM for randomly chosen states 
are collected in Table [T] As explained in the table caption, we report data for four samples of random states: 
an unbiased set of mixed states, a biased high-purity set of mixed states, a set of pure states, and a sample 
of maximally entangled states. 

The standard deviation for each quantity, as listed in Table[TJ is a combination of two types of errors: one 
that comes from performing a number of experiments on a fixed state and another incurred as we sample the 
random states. To identify the two different errors, one would need to perform an extremely large number 
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Table 1: Average values for three samples of 1000 randomly chosen states each. The top left table is for unbiased full-rank 
mixed states (r = 4 and M = in Sec. 0; the top right table is for biased full-rank mixed states (r = 4 and M ^ 0) whose 
purity tr|p 2 | exceeds 0.8 with high probability. The bottom left table is for pure states (r = 1); the bottom right table is for 
maximally entangled states. In each table, the first and second rows report the average number of qubit pairs that need to be 
measured before the benchmark value of D t hr = 0.1 is reached, for raw-data state estimation (RD) and maximum-likelihood 
estimation (ML) as well as for the product POM and for the SIC POM. The third rows show the performance factors r\ of d 3 2 1 ) 
that compares the efficiency of the two POMs. For the more relevant ML estimation, there is no significant advantage of the 
SIC POM over the product POM, except perhaps for the sample of unbiased mixed states which is dominated by low-purity 
states [l5| . 



Unbiased mixed states 




Biased mixed states 




RD ML 




RD ML 


prod 


1787 ± 102 1534 ±156 


prod 


1739 ±114 712 ±99 


SIC 


1349 ±67 1165 ±101 


SIC 


1293 ± 67 639 ± 73 


V 


1.33 ± 0.09 1.32 ±0.11 


V 


1.35 ±0.11 1.12 ±0.16 




Pure states 


Maximally entangled states 




RD ML 




RD ML 


prod 


1738 ± 118 472 ±73 


prod 


1728 ± 103 552 ±89 


SIC 


1284 ± 70 455 ±108 


SIC 


1276 ± 70 442 ± 72 


V 


1.35 ±0.16 1.09 ±0.29 


n 


1.36 ±0.11 1.27 ±0.29 



of simulated experiments, in excess of 10, 000, on each state and so reduce the former error to orders of 
magnitude smaller than the latter error. 

One lesson of the numbers in Table Q] is that the RD estimate is not only sometimes unacceptable, but 
it also performs worse than the ML estimate. Clearly, one should not rely at all on the RD estimate. 

The second lesson is that there is no significant difference between the efficiency of the SIC POM and the 
product POM for states of high purity, and there is not much of an advantage for the SIC POM when the 
states have low purity. The averages over samples of randomly chosen states do not confirm the observation 
made for Bell states in Sec. 

8. Conclusion 

We have carried out a numerical study that compares the performance of the product POM and the 
SIC POM in qubit-pair tomography. All evidence indicates that the product POM performs almost as well 
as the SIC POM; there is not much of an advantage in employing the SIC POM. Since it is much more 
challenging to realize the SIC POM in a laboratory experiment than the product POM, we conclude that 
the implementation of the SIC POM for tomographic purposes is hardly worth the trouble. 
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